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Abstract 

Three particles with large scattering length display a universal spectrum of three-body bound 
states called "Efimov trimers". We calculate the modification of the Efimov trimers of three 
identical bosons in a finite cubic box and compute the dependence of their energies on the box 
size using effective field theory. Previous calculations for positive scattering length that were 
perturbative in the finite volume energy shift are extended to arbitrarily large shifts and negative 
scattering lengths. The renormalization of the effective field theory in the finite volume is explicitly 
verified. We investigate the effects of partial wave mixing and study the behavior of shallow trimers 
near the dimer energy. Moreover, we provide numerical evidence for universal scaling of the finite 
volume corrections. 
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I. INTRODUCTION 



Few-body systems with resonant interactions characterized by a large scattering length 
a show interesting universal properties. If a is positive, two particles of mass m form a 
shallow dimer with energy E2 ~ —h^/{ma^), independent of the mechanism responsible 
for the large scattering length. Examples for such shallow dimer states are the deuteron in 
nuclear physics, the ^He dimer in atomic physics, and possibly the new charmonium state 
X(3872) in particlephysics ll]. In the three-body system, the universal properties include 
the Efimov effect |3|. If at least two of the three pairs of particles have a large scattering 
length \a\ compared to the range ro of their interaction, there is a sequence of three-body 
bound states whose energies are spaced geometrically between —h^/mrQ and —fi} /ma?'. In 
the limit 1/a — )■ 0, there are infinitely many bound states with an accumulation point at the 
three-body scattering threshold. These Efimov states or trimers have a geometric spectrum 
0: 

eJ"^ = -(e-2-/«°)"-"*n2/s:2/m, (1) 

which is specified by the binding momentum of the Efimov trimer labeled by n^,. This 
spectrum is a consequence of a discrete scaling symmetry with discrete scaling factor 6^^^°. 
In the case of identical bosons, sq ~ 1.00624 and the discrete scaling factor is e^^^^ ~ 22.7. 
The discrete scale invariance persists if a is large but finite, but in this case it connects states 
corresponding to different values of the scattering length. The scaling symmetry becomes 
also manifest in the log-periodic dependence of scattering observables on the scattering 
length a The consequences of discrete scale invariance and "Efimov physics" can be 
calculated in an effective field theory for short-range interactions, where the Efimov effect 
appears as a consequence of a renormalization group limit cycle |5|. 

While the Efimov effect was established theoretically already in 1970, first experimental 
evidence for an Efimov trimer in ultracold Cs atoms was provided only recently by its 
signature in the three-body recombination rate j6|. It could be unravelled by varying the 
scattering length a over several orders of magnitude using a Feshbach resonance. Since 
this pioneering experiment, there was significant experimental progress in observing Efimov 
physics in ultracold quantum gases. More recently, evidence for Efimov trimers was also 
obtained in atom-dimer scattering 0] and in three-body recombination in a balanced mixture 
of atoms in three different hyperfine states of ^Li [s], loj, in a mixture of Potassium and 
Rubidium atoms |lO| . and in an ultracold gas of ^Li atoms ll|. In another experiment with 
Potassium atoms |l2|, two bound trimers were observed. 

The observation of Efimov physics in nuclear and particle physics systems is complicated 
by the inability to vary the scattering length and one has to focus on the detection of 
excited states. While two- neutron halo nuclei could be bound due to the Efimov effect, the 
analysis of known halo nuclei does not show much promise for an unambiguous identification 
(See Ref. [l3| and references therein). Another opportunity to observe Efimov physics is 
given by lattice QCD simulations of three- nucleon systems [13] • A number of studies of 
the quark-mass dependence of the chiral nucleon-nucleon (NN) interaction found that the 
inverse scattering lengths in the relevant ^Si-^Di and ^Sq channels may both vanish if one 



extrapolates away from the physical values to slightly larger quark masses [ISHlTJ. This 
implies that QCD is close to the critical trajectory for an infrared renormalization group 
limit cycle in the three-nucleon sector. It was conjectured that QCD could be tuned to lie 
precisely on the critical trajectory by tuning the up and down quark masses separately (isj . 
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As a consequence, the triton would display the Efimov effect. More refined studies of the 
signature of Efimov physics in this case followed 19, 20l|. However, a proof of this conjecture 
can only be given by an observation of this effect in a lattice QCD simulation The first 
full lattice QCD calculation of nucleon-nucleon scattering was reported in [2l| but statistical 
noise presented a serious challenge. A promising recent high-statistics study of three-baryon 
systems presented also initial results for a system with the quantum numbers of the triton 
such that lattice QCD calculations of three-nucleon systems are now within sight [22|. For 
a review of these activities, see Ref. 



23 



Since lattice simulations are carried out in a cubic 
box, it is important to understand the properties of Efimov states in the box. The first step 
towards this goal is to understand these modifications for a system of three identical bosons. 

The corresponding modifications of the Efimov spectrum can be calculated in effective 
field theory (EFT) since the finite volume modifies the infrared properties of the system. The 
properties of two-body systems with large scattering length in a cubic box were calculated 
in Ref. Some properties of three-body systems in a finite volume have also been 

studied previously. For repulsive and weakly attractive interactions without bound states. 
Tan has determined the volume dependence of the ground state energy of three bosons up 

13, 



to 0{{a/Ly) [25|. InRefs. 



this result was extended for general systems of N bosons. 



It was used to analyze recent results for three and more boson systems from lattice QCD 
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For the unitary limit of infinite scattering length, some studies have been carried 
out as well. The properties of three spin-1/2 fermions in a box were investigated in (3l| . 
However, this system has no three-body bound states in the infinite volume. The volume 
dependence of energy levels in a finite volume can also be used to extract scattering phase 
shifts and resonance properties from lattice calculations 32|, |33|. For a recent application of 
this idea to the A(1232) resonance, see Refs. |3J, [35 . 

In a previous letter, we have investigated the finite volume corrections for a three-boson 
system with large but finite scattering length 36j. We have studied the modification of 



the bound state spectrum for positive scattering length using an expansion valid for small 
finite volume shifts and explicitly verified the renormalization in the finite box within this 
approximation. In the current paper, we present an extension of this work that avoids the 
expansion in the energy shift and can be applied to arbitrarily large shifts. Moreover, we 
extend our previous studies to the case of negative scattering lengths. Our results indicate 
that the finite volume corrections are subject to a universal scaling relation. We also provide 
a more detailed discussion of the technical details and our numerical method. 



II. THEORETICAL FRAMEWORK 

A. Lagrangian and Boson-Diboson Amplitude 

Three identical bosons interacting via short-range forces can be described by the La- 
grangian [H, [sj 

C = ^P^(tdt + -V^^ + ^Sd - ^ (cit^2 ^ ^ ^ _ ^rftrf^t^ + . . . ^ (2) 
\ 2 / 4 4 36 

where the dots indicate higher order terms. The degrees of freedom in the Lagrangian are 
a boson field ip and a non-dynamical auxiliary field d representing the composite of two 
bosons. Units have been chosen such that ^ = m = 1. This Lagrangian corresponds to 
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FIG. 1: Integral/sum equation for the boson-diboson amplitude. The narrow hnes denote the 
boson propagator and the thick hnes denote the fuh diboson propagator. 

the zero range limit with tq = 0. Our results will therefore be applicable to physical states 
whose size is large compared to tq. Corrections from finite range can be included via higher 
order terms in the Lagrangian but will not be considered here. The coupling constants g2 
and will be determined by matching calculated observables to a two-body datum and a 
three-body datum, respectively. The central quantity in the three-body sector of this EFT 
is the boson-diboson scattering amplitude. It determines all three-body observables and is 
given as the solution of an inhomogeneous integral equation, depicted diagrammatically in 
Fig. [TJ A detailed discussion of the infinite volume case can be found in Here we follow 



the strategy of Refs. |24l . |36| and focus on the finite volume case. 



The system of three bosons is assumed to be contained in a cubic box with side length L 
and periodic boundary conditions. This leads to quantized momenta k = ^n, n G Z^. In 
particular, the possible loop momenta are quantized. As a consequence, the integration over 
loop momenta is replaced by an infinite sum. The divergent loop sums are regulated by a 
momentum cutoff A similar to the infinite volume case. 

The finite volume modifies the infrared physics of the system but does not change the 
ultraviolet behavior of the amplitudes. Therefore, the renormalization is the same in the in- 
finite and finite volume cases. Of course, this statement only holds if the momentum cutoff A 
is large compared to the momentum scale set by the size of the volume, namely 2n/L, such 
that the infrared and the ultraviolet regime are well separated. In the numerical calculations 
presented in this work, the consistent renormalization of the results will always be verified 
explicitly. The two-body sector can be completely renormalized by matching the two-body 
coupling constant g2 to a low-energy two-body observable, namely the two-body scatter- 
ing length a or the dimer binding energy. Because of the discrete scaling symmetry, the 
three-body coupling approaches an ultraviolet limit cycle. For convenience, the three-body 
coupling constant is expressed in the form g^ = —QglH{K)/h?. The cutoff dependence of 
the dimensionless function H{A) is then given by [l| 

^^^^ ^ cos [sq log(A/A^) + arctangp] 
cos [solog(A/A*) — arctanso] ' 

where sq ~ 1.00624 . . . and A* ^ 2.62^;* is a three-body parameter that can be fixed from a 
trimer binding energy or another three-body datum. 

For the diboson lines in the three-body equation depicted in Fig. [H the full, interacting 
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FIG. 2: The full diboson propagator (thick solid line) is obtained by dressing the bare diboson 
propagator (double lines) with bosonic loops (narrow solid lines). 



diboson propagator D has to be used. This quantity corresponds to the exact two-body 
scattering amplitude. It is obtained by dressing the bare diboson propagator, given by the 
constant 4:i/g2, with bosonic loops as shown in Fig. [21 This leads to an infinite sum that 
can be evaluated analytically. For a diboson with energy E, the propagator is 



D{E) 



327r 
If 



1 1 1 



L ^ \3\ 



(4) 



In the limit L — )■ oo, this expression reduces to the full diboson propagator in the infinite 
volume case. 

Having obtained the full diboson propagator, the integral equation for the boson-diboson 
scattering amplitude can be written down explicitly. Using the diagrammatical representa- 
tion in Fig. [Hand the Feynman rules derived from the effective Lagrangian ([2]), we obtain: 



A{vME) 



9l 
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(5) 
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D{E-qo + qy4) 
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D{E-qo + q'/4) 



A{q,k;E) 



Here, p {k) are the momenta of the incoming (outgoing) dibosons, while the momenta of 
the incoming (outgoing) bosons are —p {—k). The boson legs have been put on shell but 
the diboson legs remain off-shell. The total energy of the system, E, can be treated as a 
parameter of the equation. The integration over the loop energy can be performed by virtue 
of the residue theorem. This yields the sum equation 



with 



A{p,k;E) = -^ZE{p,k) + ^ ZE{p,q)rE{q)A{q,k;E), 



ZE{p,k) 

Mq) 



p^ + p ■ k + k^ 



E 



HjA) 

A2 



3g2 



- E + V yj-e-l^^ V 
^ L\j\ 
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(6) 



(7) 
(8) 
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If the energy E is near a trimer energy E^""^ < 0, the amphtude A exhibits a simple pole 



and the dependence on the momenta seperates: 



A{p,k-E) 



E-E. 



(n) 



as E E. 



in) 



(9) 



Matching the residues on both sides of Eq. iQ, the bound-state equation 



(10) 



is obtained. Values of the energy E, for which this homogeneous sum equation has a solution, 
are identified with the energies eI^^^ of the trimer states. 

B. Cubic symmetry 



In the infinite volume case, only s-wave bound states are formed. However, in a finite cubic 
volume, the spherical symmetry of the infinite volume is broken to a cubic symmetry. In 
the language of group theory, the infinitely many irreducible representations of the spherical 
symmetry group 0(3) are mapped onto the five irreducible representations of the cubic 
group O. The representations of the spherical symmetry group are now reducible and can 
therefore be decomposed in terms of the five irreducible representations of the cubic group. 
On the other hand, a quantity ipg transforming according to the irreducible representation s 
of O can be written in terms of the basis functions of the spherical symmetry, i.e. the 
spherical harmonics Yim, via 



(11) 



where f = f/\r\ and t is an additional index needed if the representation labeled by i appears 
in the irreducible representation s more than once. The linear combinations Kgu of spherical 
harmonics are called "kubic harmonics" HT 
for values of 



as large as 12 [38 



The values of the coefficients C^^)^ are known 



In order to make contact with the infinite volume formalism, Eq. (|T0|) is rewritten using 
Poisson's resummation formula. This identity states that an infinite sum over integer vectors 
may be replaced by an infinite sum over integer vectors of the Fourier transform of the 
addend, i.e. 



E /(^) = E /(^) 



(12) 



neZ3 



meZ3 



where f{rn) = J^^ d^y e^'^'^'^'^ f{y) is the Fourier transform of /. Applying the identity ([12 
to Eq. ([TOD yields 



(13) 
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The term with n = gives the corresponding equation in the infinite volume, while the other 
terms of the infinite sum may be viewed as corrections due to momentum quantization and 
the breakdown of the spherical symmetry. The explicit recovery of the infinite volume term 
is useful for bound states, since in this case the analytic structure of the amplitude, i.e. the 
pole at the binding energy, is identical and only the position of the pole is changed. This 
approach may be inappropriate for states belonging to the continuous scattering region of 
the infinite volume. 

In order to evaluate the angular integration in Eq. f[T^ . all quantities with angular 
dependence are expanded in terms of the basis functions of the irreducible representations 
of 0(3), namely in spherical harmonics. The amplitude J-" itself is assumed to transform 
under the trivial representation Ai of the cubic group, since the i = representation is 
solely contained in Ai. The amplitude can therefore be written as 

(Ai) e 

e,t e,t m=-e 

The sum runs over those values of i associated with the Ai-representation of O. The first 
few values are £ = 0, 4, 6, 8, ... . Since the first occasion where an £-value appears more than 
once is i = 12, the summation over the multiplicity and hence the index t will be suppressed 
in the following. 

The only angular dependence of the quantity ZE{p,y) is on the cosine of the angle 
between p and y. Therefore, ZE{p,y) can be expanded in Legendre polynomials Pi^cosOpg). 
These polynomials can in turn be expressed in spherical harmonics via the addition theorem, 
yielding 

oo oo . £ 

ZE{p,y) = Y,z^E{p.y)P^{^ose.^y) = Y.^E\p.y)i^^ E ^LmUy)- (is) 

1=0 1=0 m=-£ 

The exponential function in Eq. f|T3l) can be rewritten using the identity 

oo e 

^^Ln.y ^ 4^ ^ i'UL\n\y) E Y,Uh)YeUy), (16) 

e=o m=-l 

where is the spherical Bessel function of order i. 

With these expansions, the angular integration in Eq. (fT3l) can be performed. The integral 
over the three spherical harmonics depending on y yields Wigner 3-j symbols. Projecting 
on the £th partial wave results in an infinite set of coupled integral equations for the quan- 
tities Ff-. 
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|n| = 


vji \ {Yem + Ye,' 

-n 




n 


^ = 0,m = 





i = A,m = 


£ = 4, m = 4 


101 


606 






96.23 


81.33 


1001 


9126 






-94.05 


20.69 


3501 


39678 






118.67 


-15.73 



TABLE I: Integer vector sums over real combinations of for a given absolute value n ("angular 
sums"). For higher partial waves, the value of the sum is smaller and the contribution of the partial 
wave is suppressed. 



The i' sum runs over the partial waves associated with the Ai representation. The quantity 
z'e' {Pi y) '^^'^ be calculated from Eq. ( lT5l) to be 



zg)(p,y) = (2£+i: 



py \ py J 



(18) 



Here, Qe is a Legendre function of the second kind. The three-body contact interaction 
contributes only to the s-wave, as expected. 

Since the bound states in the infinite volume are s-wave states, Eq. f|T7j) is specialized to 
the case i = 0, yielding 
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fii^I? £'=4,6,... m'=-e 



The specialization of Eq. f lTSj) to the case £ = reads 



zf{p,y) 



1 . f p^ + py + y"^ 



2py 



log 



E 



p2 _ py y2 _ 



A2 ■ 



(20) 



The second line of Eq. ( lT9l) indicates admixtures from higher partial waves. Since the leading 
term in the expansion of the spherical Bessel function is l/(L|n|?/), these contributions are 
suppressed by at least a/L. They will therefore be small for volumes not too small compared 
to the size of the bound state. Moreover, contributions from higher partial waves will be 
suppressed kinematically for shallow states with small binding momentum. This is ensured 
by the spherical harmonic in the second line of Eq. (fT9|) . For higher partial waves, the angular 
sum yields smaller prefactors relative to the s-wave {£ = 0) for a given absolute value n (see 
Table H]). Only for small lattices, i.e. when a/L is large, this behavior is counteracted by 
terms stemming from the spherical Bessel function and higher partial waves may contribute 
significantly. The first numerical calculations presented in this work were therefore mainly 
performed neglecting contributions from higher partial waves. For two specific examples, 
however, we explicitly calculate the contribution from the next partial wave i = 4 and 
demonstrate that it is small. Some details on the numerical solution of Eq. f lT9|) are given 
in the Appendix. In the following section, we present our results. 
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III. RESULTS AND DISCUSSION 



By employing the formalism laid out in the previous section, we have calculated energy 
levels in finite cubic volumes of varying side lengths. The results of these calculations are 
presented in the following. For convenience, the dependence of the energies on the boson 
mass m is reinstated in this section. 

A. Positive scattering length 

We first present results for systems with a > 0. In this regime, a physical diboson state 
with a binding energy = — l/(ma^) exists. This energy is therefore identical with the 
threshold for the break-up of a trimer into a diboson and a single boson. We choose states 
with different energies in the infinite volume, including shallow as well as deeply bound 
states: 



la: 


Er 


= -1.18907/(ma2), 


A^,a = 


5.66 


lb: 


ET 


= -27.4427/(ma2), 


A*a = 


5.66 


Ic: 


ET 


= -9440.91/(ma2), 


A^,a = 


5.66 


II: 


ET 


= -5.04626/(ma2), 


A*a = 


1.66 


III: 


ET 


= -11.1322/(ma2), 


A^,a = 


3.66 



Here, is the trimer energy in the infinite volume. Note that the states la, lb and Ic 
appear in the same physical system chracterized by A^,a = 5.66. For each of these states, 
its energy in the finite cubic volume has been calculated for various values of the box side 
length L. In order to check the consistency of our results, the calculation was carried out 
for several cutoff momenta A. For each cutoff, the three-body interaction parameterized 
by H{A) has been adjusted such that the infinite volume binding energies are identical 
for all considered values of A. If our results are properly renormalized, the results for the 
different cutoffs should agree with each other up to an uncertainty of order 1 / ( Aa) stemming 
from the finiteness of the cutoff. 

The results for the states II and III are depicted in Fig. [3] for box sizes between L = 6a 
and L = a/2. The values obtained for different cutoffs indeed agree with each other within 
the depicted uncertainty bands. Note that these bands do not represent corrections from 
higher orders of the EFT. For both states, the infinite volume limit is smoothly approached. 
As the volume becomes smaller, the energy of the states is more and more diminished. This 
corresponds to an increased binding with decreasing box size. 

In the infinite volume, state III is more deeply bound than state II. Naively, one therefore 
expects the former to have a smaller spatial extent than the latter. The size of the state can 
be estimated via the formula (— mi?^)~^/^, yielding 0.45a for state II and 0.3a for state III. 
Hence, a given finite volume should affect state II more strongly than the smaller state HI. 
This behavior can indeed be observed. For example, considering a cubic volume with side 
length L = 2a, the energy of state II deviates 10% from the infinite volume value, whereas 
the corresponding difference for state HI is less than one percent. On the other hand, the 
box size for which the energy shift of the state HI amounts to 10% is roughly 1.3a. 
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FIG. 3: Variation of the trimer energy E3 with the side length L of the cubic volume for the states 
II (left) and III (right). Plotted are three datasets for different values of the cutoff parameter A, 
together with the l/(Aa) bands. The point a/L = corresponds to the infinite volume limit. 
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FIG. 4: Variation of the trimer energy £"3 with the side length L of the cubic volume for state II. 
Plotted are three datasets for different values of the cutoff parameter A, together with the l/(Aa) 
bands, and one dataset obtained using a Taylor expanded version of the integral kernel. 



For both states, we can now form the dimensionless number r = —mE^LlQy^, where Liq% 
is the box size at which the energy differs by 10% from the infinite volume value. This yields 
r ^ 20 for state II and r ^ 19 for state III. The approximate equality of the two values of 
r may indicate the presence of universal scaling in the finite volume version of the Effective 
Theory. 

In Fig. m the three datasets obtained for state II are plotted again, this time in coniparison 



to the results from a calculation using the expansion of the kernel used in Ref. [36[ and 
described in Appendix |X1 For large values of the box size L, the results of both calculations 
agree with each other. For volumes smaller than 1.7a, the result of the calculation with 
the expansion deviates from the result of the full calculation. For this size of the volume, 
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FIG. 5: Variation of the trimer energy with the side length L of the cubic volume for the 
states la (upper curve) and lb (lower curve). Plotted are three datasets for different values of the 
cutoff parameter A, together with the l/(Aa) bands (circles: A = 200 a~^; squares: A = 300 a^^; 
diamonds: A = 400 a^^). The solid line depicts the diboson energy. 



the full result differs by about 20% from the infinite volume binding energy. Accordingly, 
the expansion employed for the integral kernel can not be applicable any longer. In Fig. SI 
results are only depicted for volumes with L > a. For smaller volumes, the results of 
calculations using the expansion become cutoff dependent, and are hence not are not properly 
renormalized anymore. 

In Fig. O we show our results for the two states la and lb in the same physical system 
characterized by A^,a = 5.66. The volume dependence of the two states is again shown for 
box sizes between L = 6a and L = a/2. The curve corresponding to the more deeply bound 
state lb shows a behavior similar to the one observed for the states II and III. The binding 
energy remains constant until the volume is small enough to affect the state. At this point, 
the energy of the state is more and more diminished as the volume becomes smaller. 

The behavior of the shallow state la is different. In the region L a, the binding is 
not further increased. The results for smaller volumes show a sharp rise of the three-body 
energy. The energy of the state becomes positive near L = 0.67a. State la is close to the 
threshold for boson-diboson scattering in the infinite volume located at Ed = — l/fma^). 



For comparison, we calculated the energy of the physical diboson according to [2J]. The 
resulting curve is the solid line in Fig. |5l Like the energy of the three-body bound states, 
the energy of the diboson is diminished in finite volumes. For volumes of the size L ^ 1.2a, 
the energy of the three-body state la becomes larger than the diboson energy and starts 
to grow. This behavior is consistent with the observation that states are always shifted 
away from the threshold in a finite volume. In the two-body sector, for example, continuum 
states have been shown to have a power law dependence on the volume, while the volume 
dependence of bound states is dominated by exponentials jl^]. The data shown for state la 
can be explained by an exponential for L > a, which characterizes the state as a bound state. 
For L < a, the data is consistent with a power law, indicating the state indeed behaves like 
a boson-diboson scattering state if its energy is above the diboson energy. 

The other investigated states do not show such a transition since their energy is well 
below the diboson energy for all considered volumes. It is unclear whether other states 
would show a similar behavior for smaller box sizes. If this is not the case, the transition 
from bound to unbound would occur only for states with infinite volume binding energies 
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1 


200 


-4.57097 


-4.57753 


-5.07581 


300 


-4.60298 


-4.60872 


400 


-4.59579 


-4.60056 


0.9 


200 


-4.31223 


-4.27927 


-6.068 



TABLE II: Energies E3{L) of state la for box sizes near L = a calculated with and without 
admixture of the £ = 4 amplitude. The dimer energy E2 is shown for comparison. 



L/a 


Aa 


E^{L)ma^, £ = 0,4 


Es{L)ma'^, i = 


1 


200 


-11.86 


-11.15 


400 


-11.79 


-11.08 


0.7 


200 


-19.06 


-20.70 


400 


-18.97 


-20.64 



TABLE III: Energies E^[L) of state II for box sizes L = a and L = 0.7a calculated with and 
without admixture of the £ = 4 amplitude. 

up to a critical value. For a conclusive analysis of the nature of the described transition, 
more data is still needed. It would be interesting to see whether such a transition shows up 
in lattice data for a state that is very close to the diboson threshold in the infinite volume. 

Since the size of the finite volumes where the state crosses the diboson energy is compa- 
rable to the size of the state itself, the breaking of the spherical symmetry may be a relevant 
effect here. To assess the influence of the higher partial waves, we extract from Eq. (fT7|) 
two coupled equations for the s-wave amplitude Fq {£ = 0) and the i = 4 amplitude F4. 
These coupled equations are then solved in a coupled channel approach. The energies of 
state la obtained by this method are summarized and compared to the s-wave only results 
in Table [ITl For L = 1.25a, the state is still below the dimer state. The inclusion of the 
higher partial wave leads to a small downward shift in the energy. For L = a, we have done 
calculations using 3 different cutoffs. The results for different cutoffs agree to two significant 
digits indicating the results are properly renormalized, but the binding is slightly reduced 
by the i = 4 contribution. For L = 0.9a, the effect of the higher partial wave is again a 
small downward shift. All results show only a deviation of about 1% from the s-wave only 
result. In summary, we find that the correction from the i = 4 admixture is extremely small. 
Moreover, the correction is of the same order of magnitude as the finite cutoff corrections 
in our calculation and a more quantitative statement requires improved numerical methods. 
In order to establish an estimate of typical corrections from higher partial waves, we have 
also performed calculations including the i = 4 contributions for the more deeply-bound 
state II. The resulting numbers are given in Tab. IIIII The investigated box sizes are about 
three times larger than the state itself. The contribution of the higher partial wave is now 
several percent. This is still a small correction but considerably larger than the finite cutoff 
uncertainty. This suggests that the extremely small corrections for state la are related to 
its unusual behavior. The dominance of the s-wave my be associated with the closeness of 
the state to the threshold. A more detailed investigation of higher partial waves including 
a description of the numerical methods will be the subject of a future publication. 
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Aa = 4500 


\n — 
IVCL — 


Ann 


J, In 

IJ (Jj 




^rel 




Orel 


OO 


-9440.91 




-9401.32 




1 


-9440.91 


0% 


-9401.36 


^ 10-6 


0.75 


-9440.91 


0% 


-9400.53 


10"^ 


0.5 


-9440.91 


0% 


-9399.38 


0.02% 



TABLE IV: Energy of state Ic for different volume sizes L and two different cutoffs A. The 
relative deviation from the infinite volume value drei = {E^ — E^)/E^ is also given. 

Since the shallower state la is more affected by a given finite volume than the deeper 
state lb, the ratio of the energies of the states is changing. In the infinite volume, this 
ratio is 23.08. For L = 1.5a, just before the shallow state crosses the dimer energy, this 
ratio has decreased to 7.4. Note that this ratio differs from the discrete scaling factor 
exp(27r/so) ~ 515 even in the infinite volume limit. This behavior is expected for shallow 
states close to the bound state threshold [if. The ratio exp(27r/so) will be approached when 
deeper states are considered. For example, the infinite volume ratio of state lb and the much 
more deeply bound state Ic is 344 and already closer to the discrete scaling factor 515. 

If we assume that the combination r = —mE^L^Qty^ is indeed a universal number, we 
are able to predict Liq% for the states la and lb. The results are Liq% ^ 0.85a for state lb 
and -^10% ~ 4a for state la. The energies calculated for these volumes are E^{L = 0.85a) ~ 
29.1/(ma^), corresponding to an 8% shift, for state lb and E^lL = 4a) ^ 1.285/(ma^), also 
corresponding to an 8% shift. These findings support the assumption that the finite volume 
corrections obey universal scaling relations. 

The states la and lb appear in the same physical system characterized by A*a = 5.66. 
In this system, an even more deeply bound state, denoted as state Ic, is present. In the 
infinite volume, the energy of this state is = — 9401.32/(ma^). This corresponds to a 
binding momentum of 97a~^. Since this is already comparable to the momentum cutoffs of 
a few hundred inverse scattering lengths employed before, we also performed calculations 
for this state using a much larger cutoff of 4500 a~^. The three-body force for this cutoff has 
been fixed such that the energy of the most shallow state la is reproduced. The resulting 
energy of state Ic is then — 9440.91/(ma^). This differs from the energy obtained using 
the smaller cutoff by 0.4%. This difference can be attributed to effects stemming from the 
finiteness of the cutoff. We have calculated the energy of this state in finite volumes of the 
sizes L = a, L = 0.75a and L = 0.5a. The results of this calculations are summarized in 
Table IIVI The values obtained using the large cutoff Aa = 4500 show no effect of the finite 
volume at all. From the infinite volume binding energy, the size of the state can be estimated 
via (mi?|°)~^/^ to be O.Ola. So, we do not expect any visible effect since the finite volume 
is fifty times as large as the state itself. However, for the smaller cutoff Aa = 400, there 
are very small deviations from the infinite volume energy. But these deviations are smaller 
than the uncertainty stemming from the finiteness of the cutoff, which is estimated to be of 
order l/(Aa) = 0.25%. 
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a/L a/L 

FIG. 6: Variation of the trimer energy with the side length L of the cubic volume for the states 
NI (left) and Nil (right). Plotted are three datasets for different values of the cutoff parameter A, 
together with the l/(Aa) bands. The point a/L = corresponds to the infinite volume limit. 



B. Negative scattering length 

Now we turn to systems with a < 0. In this regime, the two-body interaction is attractive 
but no diboson bound state exists in the infinite volume limit. The only possible breakup 
process for a three-boson bound state is therefore the breakup into three single bosons. The 
threshold for this process is E = 0. As for the case with positive scattering length, we choose 
states with different energies in the infinite volume: 

NI: = -9/(ma2), A,a = 5.55, 

Nil: E^ = -4:/{ma'^), A,a = 4.23, 

NIII: E^ = -0.2/(ma2), A,a = 2.60, 

Here, E^ is the trimer energy in the infinite volume. For each of these states, its energy 
in the finite cubic volume has been calculated for various values of the box side length L. 
In order to check the consistency of our results, the calculation was carried out for several 
cutoff momenta A. For each cutoff, the three-body interaction parameterized by H{A) has 
been adjusted such that the infinite volume binding energies are identical for all considered 
values of A. 

The results for the states NI and Nil are shown in Fig. [H] for box sizes between L = 2a 
and L = a/2. The results for the state NIII with box sizes between L = 8a and L = a/2 
are depicted in Fig. [71 The values obtained with different cutoffs all agree with each other 
within the uncertainty bands indicating proper renormalization. Our findings are similar 
to those in the positive scattering length regime. All three states smoothly approach the 
infinite volume limit. As the box size becomes smaller the energy of the state is more and 
more diminished. The overall behavior is identical to the one in the positive scattering 
length case described in the previous section. 

The more deeply bound a state is in the infinite volume, the smaller is its spatial extent. 
Estimating the size via the formula (— mi?|°)~^/^ yields a/3 for state NI, a/2 for state Nil, 
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FIG. 7: Variation of the trimer energy E3 with the side length L of the cubic volume for the states 
NIII. Plotted are three datasets for different values of the cutoff parameter A, together with the 
l/(Aa) bands. 



State 




£'3(L = a) ma? 




^10% /a 


^100% /a 


NI 


-9 


-11.6 


29% 


1.2 


0.7 


Nil 


-4 


-8.7 


118% 


1.8 


1.05 


NIII 


-0.2 


-6.7 


3250% 


7.2 


4.2 



TABLE V: For the states NI, Nil and NIII, the infinite volume energy and the energy in 
a finite volume with side length L = a are shown together with the relative deviation from the 
infinite volume value d^f^i = {E3 — E^) / E^ . Also given are the box sizes leading to an energy shift 
of 10% and 100%, respectively. 



and 2.2a for state NIII. A given finite volume should therefore affect state Nil more than 
state NI, and state NIII should be the most affected. When considering a cubic volume with 
side length L = a, we find the energies given in Table El The relative deviation of state NI 
is indeed four times smaller than the shift for state Nil and a hundred times smaller than 
the shift for the shallow state NIII. On the other hand, the box length Liq%, for which the 
energy of each states deviates 10% from its infinite volume value, is smaller the more deeply 
bound a state is. These values are also given in Table |Vl 

As in the case of positive scattering length, we have investigated the dimensionless com- 
bination r = —mE^L^QO/^. From state NI, we obtain r ^ 13. This leads to the predic- 
tions LiQ% = 1.8a for state Nil and Lio% = 8a for state NIII in good agreement with 
the explicitly calculated values. Additionally, we form another dimensionless combination, 
r' = —mE^ Ll^Qty^, where Liqq% is the box side length where the energy of the state is twice 
as large as the infinite volume value. For the three investigated states, this box size is also 
given in Table |Vl From the value for the state NI, we get r' ^ 4.4. From this value, we 
predict -^100% = l-05a for state Nil and Liqq% = 4.7a for state NIII. Again, the value for 
state Nil is as predicted, whereas the value for state NIII is about 10% off the prediction 
from universal scaling. 
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For states far away from the threshold, we expect that the regimes of negative and positive 
scattering length are governed by the same scaling factor. Therefore, the dimensionless 
combination r should, for such deeply bound states, have a common value for both signs 
of a. For positive scattering lengths, state lb is an example of a rather deeply bound state. 
For this state, we find Liq% = 0.79a and from that r = 16.8. For negative scattering lengths, 
we chose a state with = —2'!/ {ma?) by setting A^a = 8.54. The energy of this state is 
shifted by 10% in a volume with side length 0.73a, yielding r = 14.4. The values of r for the 
two signs of the scattering length are indeed close to each other for these two states. This 
behavior provides numerical evidence for the universality of finite volume effects for both 
positive and negative scattering lengths. 



IV. CONCLUSIONS 



In this paper, we have extended our earlier studies of the Efimov spectrum in a cubic box 
with periodic boundary conditions 3^. The knowledge of the finite volume modifications 



of the spectrum is important in order to understand results from future 3-body lattice 
calculations. Using the framework of EFT, we have derived a general set of coupled sum 
equations for the Efimov spectrum in a finite volume. 

Speciahzing to £ = 0, we have calculated the spectrum for both positive and negative 
scattering lengths and verified the renormalization in the finite box explicitly. We have 
removed the expansion for small finite volume shifts used in Ref. [33] and presented an 
extension that can be applied to arbitrarily large shifts. Typically, the binding of all three- 
body states increases as the box size is reduced. Moreover, we provided a more detailed 
discussion of the technical details and our numerical method. We have investigated the 
breakdown of the linear approximation in detail and find that the expansion is applicable as 
long as the finite volume shift in the energy is not larger than 15-20% of the infinite volume 
energy. 

We have studied the spectrum for positive and negative scattering lengths in detail and 
provided numerical evidence for universal scaling of the finite volume effects. The scaling 
properties can be quantified by the dimensionless product r of the three-body binding energy 
in the infinite volume limit and the square of the box length corresponding to a finite volume 
shift of ten percent of the infinite volume energy. For sufficiently deep states, we obtained 
numerical evidence that r approaches a universal number for both signs of the scattering 
length. These findings suggest that the properties of deeper states in the finite volume can 
be obtained from a simple rescaling and do not require explicit calculations. A more detailed 
analysis of this issue along the lines of 3^ is in progress. 



For positive scattering lengths, we have investigated a spectrum of three states in the 
same physical system and studied the behavior of the shallowest state in the vicinity of the 
dimer energy which specifies the scattering threshold in the infinite volume. We found that 
this state drastically changes its behavior as a function of the box length once its energy 
becomes equal to the dimer energy. At this point the energy of the shallowest three-body 
state starts to grow and eventually becomes positive. The observed behavior is consistent 
with exponential suppression of the finite volume corrections below the dimer energy and 
power law suppression above. 

The effect of an admixture of the £ = 4 partial wave has been investigated for two different 
three-body states for box sizes of the order of the size of the state. For a generic state well 
separated from threshold, we found corrections of the order of a few percent for volumes 



16 



that are about three times larger than the state itself. For the state closest to threshold, the 
effect of this admixture turns out to be surprisingly small for volumes about twice as large 
as the state and was found to be less than 1%. As the finite cutoff corrections for this state 
are about the same size, a more quantitative study requires an improved treatment of these 
corrections. Such studies are in progress. 

Finally, our method should be extended to the three-nucleon system. This requires also 
the inclusion of higher order corrections in the EFT and finite temperature effects as lattice 
calculations will inevitably be performed at a small, but non-zero, temperature. Work in 
these directions is in progress. With high statistics lattice QCD simulations of three-baryon 
systems within reach [22|, the calculation of the structure and reactions of light nuclei 
appears now feasible in the intermediate future. Our results demonstrate that the finite 
volume corrections for such simulations are calculable and under control. This also opens 
the possibility to test the conjecture of an infrared limit cycle in QCD for quark masses 



slightly larger than the physical values [18 
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Appendix A: Numerical treatment 

In this appendix, we give some details on the numerical solution of Eqs. f ITTt [T^ . The 
starting point for a first numerical treatment of the formalism is the homogeneous integral 
equation 

FM = i rivv'z'°'{p.yMv)(i + J: ?^^)fm. (Ai) 

The first step is to transform this equation into a finite-dimensional problem. Due to the 
oscillatory nature of the integrand, it is not sensible to use a finite number of sampling points. 
Therefore, a set of basis functions will be used. The choice of the basis functions is guided 
by the knowledge of the bound-state amplitude in the infinite volume case. Asymptotically, 
the amplitude in the infinite volume behaves like [l|, [HI 

Fo{p) -> - cos (so log(p/p*)) , 
P 

with the universal number sq = 1.00624. . . and a momentum scale p^,. The bound-state 
amplitude obtained by the infinite volume formalism has precisely such a form, with only 
a few zeros in the interval [0;A] (See, e.g., Ref. jsj). Therefore, the basis functions are 
chosen as Legendre functions Pi with logarithmic arguments as 

e.(p) = -P.(21og2(p/A+l)-l). (A2) 
P 
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These basis functions are orthogonal with respect to a suitably chosen scalar product: 



dp - 

P + A 



2i + l 



-Si 



(A3) 



The amplitude Fq in Eq. (lAip is replaced by its expansion in the basis functions X]^o fj^jip) 
and the ith component is projected out using the scalar product. The resulting equation 



log 2 



2i + 



(A4) 



with 



K,,{E) 



TT 



\ P" 
dp r 

P + A 



L\n\y 



can be interpreted as a matrix equation when truncating the set of basis functions. 
A non-trivial solution of Eq. ( ]A4I) can be found if the condition 



det 



J 2i + l ^ 







(A6) 



is fulfilled. Values of E, for which this is the case, are then considered as trimer energies in 
the finite volume. This is only a valid interpretation if the truncation of the basis and the 
neglect of the higher partial waves induce only small corrections for the result. 

The value of the parameter E for which the condition flA6p is satisfied is found via a 
root finding algorithm. This requires the calculation of Kij{E) in each iteration. To save 
numerical effort, it is possible to expand the kernel around the binding energy in infinite 
volume if the shift in the binding energy is small. The expansion is done up to first order: 



K,,iE) ^ K,,iE^) + 



dE 



(E-E' 



for 



E=E° 



E-E" 



< 1. 



(A7) 



This expansion has been used to obtain the results in |36|]. It is applicable as long as the 



finite volume shift in the energy is not larger than 15-20% of the infinite volume energy. 

In the following, details on the numerical methods used to calculate the kernel matrix 
elements as defined in Eq. flASP are presented. The integration over p is performed using a 



logarithmically distributed Gauss- Legendre quadrature with 64 points. 

The integrand of the y-integration is strongly oscillating due to the sin{L\n\y) term. 
Therefore, this integral is evaluated using the Fast Fourier Transformation technique (FFT). 
The computation of Fourier-type integrals via FFT is explained in detail in ji^]. The 
sampling for the FFT is the most time-consuming part of the calculation. The number of 
points to sample, M, is determined by the highest "frequency" that is to be accessed. The 
frequencies in the present case are the values of L\n\. The highest possible frequency is 
connected to M and the integration range A by 



I maxi ^ 



M > 



-^I'^maxIA 
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(A8) 
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FIG. 8: Convergence of the sum over vectors in J? as a function of the radius R. 

For realistic values of these parameters, namely box sizes of a few a, an |nmax| of about 3000 
and cutoff values of several hundred a~^, a typical value for M is 2^'^ or about one million 
sampled points. 

A summation over three-dimensional integer vectors of a quantity depending only on the 
absolute values involves substantial double counting when naively done. To circumvent this 
problem, an ordered list with all absolute values that are possible for vectors in 1? and their 
multiplicity has been created. The result of the summation can be viewed as converged 
when going to absolute values of about R ~ 5000 (see the example in Fig. |H]). For radii R 
from about 2000 on, the intermediate sums oscillate around the converged result. To reduce 
runtime, the result of the summation has been calculated as a mean of 15 intermediate sums 
for radii from 2750 to 3500. 

The calculation of the matrix elements has been performed on the cluster of the HISKP 
at the University of Bonn. As already stated above, the number of sampling points for the 
FFT is the parameter with the largest influence on the runtime. When using the expanded 
integral kernel, the matrix and the derivative matrix have to be calculated only once. For a 
typical sample size of 2^°, this takes about 100 minutes. When going to smaller volumes, the 
shifts in the binding energy become larger and the Taylor expansion of the integral kernel 
breaks down. In this case, it is inevitable to recalculate the kernel matrix in each iteration 
of the root finding algorithm. This amounts to 10 to 15 evaluations of the kernel matrix and 
a runtime of 8 to 10 hours per data point. 
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